Upper critical field for electrons in two— dimensional lattice 
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We address a problem of the upper critical field in a lat- 
tice described by a two-dimensional tight-binding model with 
the on-site pairing. We develop a finite-system-approach 
which enables investigation of magnetic and superconducting 
properties of electrons on clusters, consisting of a few thou- 
sand sites. We discuss how the quasiparticle density of states 
changes with the applied external magnetic field and present 
the temperature dependence of the upper critical field. We 
also briefly discuss possible extension of the model to account 
for the properties of high-temperature superconductors. 

74.25.Ha, 74.60. Ec, 71.70.Di 



The issue of the critical field consists of two different 
phenomena, namely a movement of electrons in a pe- 
riodic potential under influence of a magnetic field and 
superconductivity. Each of these phenomena has been in- 
vestigated since many decades and many solutions in lim- 
iting cases are known at present. Concerning the move- 
ment of electrons one deals with two limiting cases: free 
or nearly free electrons in a magnetic field, when Lan- 
dau levels structure sets on , and electrons in a peri- 
odic potential in the absence of magnetic field, when the 
solutions are Bloch waves which lead to energy bands. 
Away from these limiting cases the situation is much 
more complicated. Application of magnetic field to the 
two-dimensional (2D) electron system in tight-binding 
approximation leads to a fractal energy spectrum known 
as Hofstadter's butterfly, where very small changes in 
magnetic field can result in a completely different spec- 
trum H H . Electrons on a lattice are gauge-invariantly 
coupled with a U(l) gauge field by introducing phase 
factors in the kinetic-energy hopping term, i.e., the wave 

function acquires a factor exp f |§ A • dlj , where A 

is the external classical vector potential, when an elec- 
tron hops from site i do site j ||. The Zeeman term is 
neglected. The same energy spectrum can be obtained 
in a nearly-free-electron method with a weak periodic 
perturbation introduced into the Landau-quantized 2D 
electron system HJzlll ■ 

On the other hand, the influence of a magnetic field on 
superconductivity is usually described by phenomenolog- 
ical Ginzburg-Landau theory Q , (or Lawrence-Doniach 
theory in the case of layered superconductors Jl(J ) where 
the magnetic field is treated semiclassically. This ap- 
proach was later justified also at the microscopic level 



M, but the temperature dependence of physical quanti- 
ties is also of a phenomenological character and therefore 
its validity is limited. 

Although, there is a general agreement that exter- 
nal magnetic field reduces the critical temperature, pos- 
itive curvature of H C 2 (T) observed in high-temperature 
superconductors is still a matter at issue. The 
most important differences between standard BCS-type 
and high-temperature superconductors are related to the 
presence of strong electronic correlations and specific ge- 
ometry of high-T c materials. In this paper we address an 
important problem concerning the upper critical field for 
electrons described by the two dimensional tight-binding 
model. 

Our starting point is two-dimensional square lattice 
immersed in a perpendicular, uniform magnetic field. 
The mean-field Hamiltonian is of the form 



Here, c\ a (c icr ) creates (annihilates) an electron with spin 
<7 on site i, V stands for the magnitude of the on-site 
attraction and A is the vector potential corresponding 
to the external magnetic field H. Similarly to Gor'kov's 
approach we introduce local superconducting order pa- 
rameter IpJ 
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A > = ( c U c *t)' 



(2) 



which, in general, can be site-dependent. According to 
the Peierls substitution || the original hopping integral 
Uj is multiplied by a phase factor, which accounts for 
coupling of electrons to the magnetic field 



hi ( A ) = hi ex P 




A-dl 



(3) 



In order to derive the self-consistent equation for the gap- 
function, we make use of unitary transformation and in- 
troduce 



= V[/t , c . 

icr / j vni Ky ' 



(4) 



where the unitary matrix U im consists of eigenvectors of 
the hermitian matrix (A) 
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This unitary transformation determines energy spectrum 
of the system in the normal state in the presence of ex- 
ternal magnetic field. In the absence of magnetic field, 
Eq. (4) represents a transformation to the momentum 
space, namely Uj m = l/Nexp(iRj ■ k m ). In the case 
of free electron gas external magnetic field leads to the 
occurrence of rotationally invariant states corresponding 
to the Landau orbits. In order to visualize the impact 
of magnetic field on electrons in the 2D lattice we have 
calculated the resulting current distribution. Within the 
framework of linear-response theory the current oper- 
ator is given by Ji(x,y) = —8H/dAi(x,y), where (x, y) 
denotes spatial coordinates and Vs are unit vectors in the 
lattice axes directions. Results obtained on a finite clus- 
ter with applied magnetic field ea 2 H z /(Tic) — 0.1 are pre- 
sented in Figure 1. Note that radiuses of corresponding 
Landau orbits are of the order of a few lattice constants 
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FIG. 1. Current distribution originating from external 
magnetic field calculated for a 20x20 cluster with fixed 
boundary conditions. Only nearest-neighbor hopping has 
been taken into account. The size of an arrow is proportional 
to the magnitude of current. 

Equations of motion lead to formally exact expression 
for the anomalous Green's function 

(w - EJ <(a mT | a nL )) = -Vj2\U? m U: n , 



(6) 



At temperatures close to T c one can restrict himself 
only to terms linear in the superconducting order pa- 
rameter. Therefore, the Green's function \( a ni^ I a n i 



which enters the above equation can be calculated in the 
normal state and one can derive a system of linearized 
gap equations which determine the upper critical field 

A =yVA-V[/ U- U* U* 

j v / j % / j .?n .?m %m m 



tanh 
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(7) 



In the absence of magnetic field one can choose the or- 
der parameter as a site-independent quantity (Aj = A). 
Then the above equation can be easily reduced to the 
standard BCS form. This formula is gauge-invariant 
and is valid for any dispersion relation determined by 
the hopping integral tjj. However, in order to simplify 
numerical calculations we restrict ourselves only to the 



nearest-neighbor hopping with ti 



t and choose the 



Landau gauge A = H z (0, x, 0). Here, H z denotes a mag- 
nitude of external magnetic field perpendicular to the 
(x, y) plane. Such a form of A neglects the effects of dia- 
magnetic screening supercurrents induced by the applied 
field what means that the true vector potential should 
be determined self-consistently. However, close to the 
transition temperature i.e., in the limit of infinite Lon- 
don penetration depth, such an approximation is correct. 
Then, the unitary matrix U in takes on the form: 

U m ee U x , v {p XlPy ) = N- 1 / 4 e lp vV a g (p x ,p y ,x) , (8) 

where (x,y) enumerate the lattice sites R J9 = e x ax + 
e y ay and p y is the wave vector in the y direction. 

It follows from Eq. (5) that g (p x ,p yl x) must fulfill 
the Harper's equation |1: 

9 (Px , P y i x + l) + 2 cos (hx - p y a) g (p x , p y , x) 

+g(p x ,p y ,x- i) = t~ l E (p x ,Py) g (p x ,p y ,x) , (9) 

where h — ea 2 H z /(hc). /i/(27r) can be interpreted as a 
ratio of the flux through a lattice cell to one flux quantum 
|| . In the absence of magnetic field p x corresponds to the 
x component of the wave vector p. As the choice of Lan- 
dau gauge breaks the translational symmetry along the x 
axis, for H ^ p x represents a quantum number which, 
however, can not be identified as a component of the wave 
vector. The choice of this gauge allows one to take Ui n in 
the form given by Eq. (8) and reduces the original two- 
dimensional eigenproblem (diagonalization of the kinetic 
part of the Hamiltonian) to a one-dimensional difference 
equation. A thorough analysis of Harper's equations can 
be found in Ref. @. 

Due to the plane-wave behavior in the y direction, 
there is a solution of Eq.(7) which does not depend ex- 
plicitly on y: 



A, = A^ = A. 



(10) 



For the specific form of Ui n , as given by Eq. (8), the up- 
per critical field is determined by the following equation: 
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^X! A * Yl 9(p x ,P y ,x')g(k x ,- Py ,x') 

x 9 {k x ,-p yl x) g (p xl p y ,x) 

tanh^^+tanh^'"^ 



2k B T 



2 '' ljT . (11) 



2[E(p x ,p y )+E(k x ,-p y )] 

In order to evaluate the upper critical field one has 
to start with solving the Harper's equation. The corre- 
sponding energy spectrum E(p x ,p y ) has been obtained 
for the first time by Hofstadter and constitutes a self- 
similar, fractal structure known as the Hofstadter but- 
terflyg. 

In the thermodynamic limit Eq. (11) is actually an in- 
finite system of equations and can not be solved exactly. 
Therefore, in order to get the first insight into the prop- 
erties of the system under consideration, we have per- 
formed numerical calculations for finite systems. In par- 
ticular, we have investigated square clusters which consist 
of a few thousand lattice sites. Due to the broken transla- 
tional symmetry, we have introduced periodic boundary 
conditions only along the y axis and fixed boundary con- 
ditions in x direction. It means that for the M x M 
cluster we have taken into account M points along x axis 
and M values of the wave vector p y G (— n/a, n/a]. 
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FIG. 2. Density of states obtained for a 50 x 50 cluster and 
h — 27r/10. Vertical bars show the energy spectrum obtained 
for infinite system (Hofstadter's butterfly). Energy levels ob- 
tained from the cluster calculations are represented by the 
Lorentz function with the width equal 0.03i. 



which occurs in cluster calculations with fixed and 
riodic boundary conditions taken in both directions 

To check the influence of finite size effects we have cal- 
culated the density of states of 2D electron lattice gas in 
the normal state. Fig. 2 shows that the density of states 
calculated for a 50 x 50 cluster reproduces very well the 
results obtained on the basis of Hofstadter's procedure 
for infinite system g. 

As a further verification of our cluster approach we 
have compared the critical temperature calculated with- 
out the external magnetic field with exact results for the 
2D lattice obtained from the BCS-type gap equation (see 
the inset in Fig. 3). The critical temperature obtained 
from the cluster calculations is always a bit lower than 
the exact value simply due to the absence of van Hove 
singularity in finite cluster. Since our method works well 
in both limiting cases i.e., in the normal state influenced 
by external magnetic field and in superconducting state 
analyzed without the field, we have used this approach 
to tackle the problem which is fundamental in the inter- 
mediate region, namely the influence of the external field 
on superconductivity. Numerical solutions of Eq. (11) 
are shown in Fig. 3. 
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FIG. 3. Reduced upper critical field obtained for 40 x 40 
and 50 x 50 clusters plotted as a function of temperature. 
V — 0.7t has been used. Filled circles and triangles denote 
results which fulfill the criterion fcgT > 8t/Al 2 . The inset 
shows the superconducting transition temperature calculated 
in the absence of magnetic field versus the magnitude of pair- 
ing potential. Here, continuous and dashed lines show exact 
result for infinite lattice and solution for 40 x 40 cluster, re- 
spectively. 



Making use of the fixed boundary conditions in the x 
direction the Harper equation (9) simplifies to an eigen- 
problem of a tridiagonal matrix with all the off-diagonal 
elements equal unity. An additional effect originating 
from such specific boundary conditions is the absence 
of unphysical degeneracy of states at the Fermi level 



However, one has to bear in mind that our cluster 
calculations are not valid in genuinely low temperatures 
when the Cooper pair susceptibility accounts only for 
very few poles of the Green's function [E (p x ,Py)] instead 
of a continuous density of states. The simplest criterion 
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of validity for a M x M cluster is an assumption that 
temperature (&bT) must be larger than an average dis- 
tance between different quasiparticle energies (~ 8t/M 2 ). 
We have found that, similarly to the free electron gas 
described by Gor'kov equations pLbJ, also in the tight- 
binding model H C 2 (T) exhibits the negative curvature. It 
is remarkable that the critical temperature is a smooth 
function of applied magnetic field despite the fact that 
the energy spectrum is strongly affected even by small 
changes of the field ||. However, for realistic values of 
magnetic field, h ~ 1CP 4 -j- 10~ 3 , the band splits into a 
huge number of subbands (oc h^ 1 ) with the gaps between 
them much smaller than fc^T. Therefore, macroscopic 
properties of the system remain smooth functions of the 
magnetic field. 

To conclude, we have investigated the relationship be- 
tween superconductivity and an external magnetic field 
in a simple model of electron gas in a 2D lattice. The pro- 
posed method allows one to analyze such a model on large 
clusters, of size of the order of a few thousand lattice sites, 
carrying out exact calculations. The main result is the 
temperature dependence of the upper critical field (or, 
equivalently, the applied magnetic field dependence of 
the critical temperature) . We have found that the exter- 
nal magnetic field suppresses superconductivity, what re- 
mains in agreement with a general filling, although there 
are papers where the reentrance of the superconductiv- 
ity due to Landau quantization in strong magnetic field 
is analyzed |l3[ ]. In contradistinction to the Gor'kov ap- 
proach, which is valid for free electron gas, we have car- 
ried out calculations for a system with different geometry 
i.e., for the 2D electron gas on a square lattice. Such a ge- 
ometry is believed to play a crucial role in high-critical- 
temperature superconductors. However, in spite of the 
differences between these two approaches the curvature 
of the transition curve H C 2{T) in both cases is negative. 
This result suggests that the lattice geometry itself is in- 
sufficient to explain the origin of the positive curvature 
of the upper critical field observed in the copper-oxide 
high T c superconductors and additional factors have to 
be taken into account. The most important feature of 
these materials neglected in our approach is the presence 
of strong correlations resulting from the on-site repulsion 
on the copper sites, which can considerably suppress pair- 
ing in the s-wave channel. On the other hand, in order to 
analyze the influence of these correlations on the upper 
critical field one has to diagonalize the full Hamiltonian 
with the interaction term, what is possible only for very 
small clusters for which the finite size effects play much 
more important role. Apart from the finite size effects, 
there are too few points in the Brillouin zones of such 
small clusters to reproduce the structure of the energy 
levels of infinite system. In the framework of the pro- 
posed approach i.e., without the necessity of the above 
mentioned diagonalization of the Hamiltonian in the full 
Hilbert space, one can take into account not only local, 
but also non-local pairing. Such an extension would lead 
to an anisotropic superconductivity, which can survive in 



spite of the presence of strong on-site repulsion |l4| . This 
problem is currently under investigation. 
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